Phase-field model of long-time glass-like relaxation in binary fluid mixtures 
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We present a new phase-field model for binary fluids, exhibiting typical signatures of soft-glassy 
behaviour, such as long-time relaxation, ageing and long-term dynamical arrest. The present model 
allows the cost of building an interface to vanish locally within the interface, while preserving 
positivity of the overall surface tension. A crucial consequence of this property, which we prove 
analytically, is the emergence of free-energy minimizing density configurations, hereafter named 
"compactons" , to denote their property of being localized to a finite-size region of space and strictly 
zero elsewhere (no-tails). Thanks to compactness, any arbitrary superposition of "compactons" still 
is a free-energy minimizer, which provides a direct link between the complexity of the free-energy 
landscape and the morphological complexity of configurational space. 

PACS numbers: 47.11.-j,47.61.Jd,47.57.jb 



The coarsening mechanisms by which a binary fluid 
mixture attains an ordered state upon a deep quench 
from a high-temperature disordered phase, continue to 
attract a great deal of attention in the scientific commu- 
nity . Apart from their paramount practical inter- 
est, these phenomena still set a fascinating challenge to 
the foundations of non-equilibrium thermodynamics, be- 
cause of their competition-driven long-time relaxation, 
often denoted as "glassiness" or " self-glassiness" . The 
literature on glassy systems is huge, covering materials 
as diverse as colloids, block copolymers, proteins, glass 
forming liquids and many others jl], HJ ■ Equally broad 
is the spectrum of theoretical/computational techniques 
employed for their study, such as replica methods, mode- 
coupling theory, mcsoscopic kinetic models, as well as 
molecular dynamics, Montecarlo and Langevin simula- 
tions More specifically, coarsening phenomena in bi- 
nary mixtures are typically described by Langevin equa- 
tions, governing the space-time evolution of the order 
parameter, i.e. the density deficit between the two fluid 
densities [J]. Depending on the specific details, different 
exponents are then predicted for the power-law growth of 
the coarsening length, the typical linear size of the coars- 
ening domains. In soft-glassy materials, however, do- 
main growth is observed to undergo long-term slow-down 
and possibly even dynamical arrest. In this Letter, we 
present a new phase-field Landau- Ginzburg (LG) model 
exhibiting most typical signatures of self-glassiness, such 
as long-time relaxation, ageing and long-term dynamical 
arrest. The present model can be analytically derived, 
bottom-up, from a mesoscopic kinetic scheme for com- 
plex fluids with competing short-range attraction and 
long-range repulsion [5|. Similarly to previous phase- 
field models p, 0|, the stiffness coefficient, controlling 
the cost of building and maintaining an interface between 
the two fluids, acquires a dependence on the local value 
of the order parameter. However, unlike any previous 
work we are aware of, instead of triggering local insta- 
bilities by sending the leading interface term to negative 



values everywhere across the interface, and then stabi- 
lizing through higher order inhomogeneities 0, here the 
stiffness becomes zero only locally within the interface, 
thereby preserving the positivity of the overall surface 
tension. This subtle difference spawns far-reaching conse- 
quences. Indeed, the present model is analytically shown 
to promote the emergence of stable, finite-support, den- 
sity configurations, which we name "compactons" . The 
dynamics of these "compactons" is then shown to be ul- 
timately responsible for the self-glassiness of the binar 
mixture. Here and throughout, at variance with ref. 
the term "compacton" is kept within quotes, to imply 
that it just refers to the property of these density excita- 
tions of being localized within a finite-support region of 
configuration space, and zero everywhere else, through- 
out the evolution. The emergence of "compactons" is 
hereby discussed analytically, both in the continuum and 
discrete versions of our phase-field models. Typical sig- 
natures of self-glassiness, such as ultra-slow relaxation, 
ageing and dynamical arrest, are further demonstrated 
by direct numerical simulations. Let us start by consid- 
ering the following LG like phase-field equation: 
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where <fi(x;t) is the order parameter, taking values (f> = 
±1 in the bulk, and (f> = at the two-fluid inter- 
face. In the above, V((f>) is the bulk free-energy den- 
sity, which we shall take in the standard double- well form 
V(4>) = —\4> 2 + i0 4 i supporting jumps between the two 
bulk phases, <j> = ±1; where r] is a white noise ^-correlated 
in space and time, with variance e. The key ingredient 
of our model rests with the specific form of the stiffness 
function D(cj)), describing the lowest order approxima- 
tion to the energy cost of building an interface between 
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FIG. 1: An example of a static gas of "compactons" for the 
case Do = and fixed E < |. The solution is constructed 
from an arbitrary superposition of stable, finite-support den- 
sity configurations, each one centered around its Xg and iden- 
tically zero outside the segment [xq < x < x ^ + l e \. The 
size l e is set by E and D2 as in eq. 

the two fluids. In the standard LG formulation, this is 
a constant parameter Do, fixing the value of the surface 
tension, through the relation 7 ~ Dq f(d x ^>) dx, x being 
the coordinate across the (flat) interface. Positive values 
of 7 promote coarsening, as a result of the surface ten- 
sion tendency to minimize the surface/volume ratio of 
the fluid. Negative values of 7, on the other hand, pro- 
mote an unstable growth of the interface, an instability 
that is usually tamed at short-scale by including higher- 
order "bending" terms of the form ~ k(A0) 2 where k 
is referred to as bending rigidity. It is readily seen that 
with Dq < and k > 0, the system undergoes instabil- 
ities, which are typically responsible for pattern forma- 
tion 0, Q . Such instabilities are then stabilized at short 
scales by a positive bending rigidity. Gompper et al., 
among others Q , studied the case with piece- wise con- 
stant D{4>) to describe microemulsions pj. Our model 
belongs to the same class as Gompper's one, with 

D{cj>) =D Q + D 2 cf> 2 (3) 

yet with a crucial twist: instead of sending Do to negative 
values, in order to trigger local interface instabilities, we 
just set D = and achieve a local zero-cost condition, 
D(4>) = 0, just at <j> = 0, by letting D 2 > 0. Thermo- 
dynamic stability of the interfaces is still secured, since 
7 > 0, and consequently we resolved to set the bending 
rigidity to k — in <j2j) , so as to single out the effect of the 
modulated stiffness D (<p) alone. In the following, we shall 
show that the peculiar feature discussed above holds the 
key for observing ultimate arrest of the fluid. As antici- 
pated, this is due to the onset of complex density config- 
urations, resulting from arbitrary superpositions of sta- 
ble, finite-support density configurations, which we name 
"compactons" . 

Let us then present our analytical analysis by looking 
at the one-dimensional, stationary solutions of the eqs. 
([HO in the limit Dq — s- and no noise (e = 0). One 



quadrature yields 

l -D 2 <p 2 (d x ct ) f+ 1 -^--^ = E (4) 

where E < 1/4 is an arbitrary constant fixing the 
energy of the configuration. Note that the relation 
between E and F for a single compacton is F — 
J (D 2 (j) 2 (d x (j)) 2 — E) dx. The analytical solution is pro- 
vided by 

</)e(x) = ±\/l - cosh(£) + e sinh (£)x(x; x , O (5) 

where £ = (x — Xo)/ld, Id = \/D 2 /2, e = 2\fE. Here 
xq is an arbitrary location and x is t ne characteristic 
function (\ = 1 inside and \ = outside) in the seg- 
ment xo < x < xo + l e , l e = Z^arctanh (2e/(l + e 2 )) (see 
figure [1]). Several comments are now in order. First, 
this solution is compact, i.e., it is identically zero outside 
the segment [xq < x < xq + l e ]. This property is cru- 
cially related to the vanishing of the prefactors in front 
of the differential operators, which allows discontinuity 
in the slope of </>(x). The location of the segment xo is 
arbitrary because of translation invariance, whereas its 
extension l e is dictated by the "energy" E. Under the 
condition that Id be real, i.e., D 2 > 0, a finite amount of 
energy E > allows the nucleation of a compacton of size 
l e > 0. The "compacton" can eventually invade the sys- 
tem, l e — > L, L being the size of the domain, a condition 
which is met at an energy value El = 1/4, since l e — > 00 
as E — > 1/4. More interesting, however, is the possibility 
of a gas of "compactons" , which can "invade" the system 
at lower values E < El, by simply superposing a collec- 
tion of disjoint compactons centered upon different values 
of xq- The possibility of such a linear superposition of 
elementary solutions of a highly non-linear field theory, 
is again a precious consequence of compactness. Since 
"compactons" do not overlap, they obey a non-linear su- 
perposition principle 4>i) n = ^ i <^>" for any power n, 
where i = 1,N labels a series of "compactons" eventu- 
ally covering the full interval, J^^Li = L. As a result, 
an arbitrary superposition of "compactons" still obeys 
the generalized LG equation. By using the above non- 
linear superposition principle, a standard stability anal- 
ysis shows that, as long as the overall surface tension is 
positive, 7 > 0, the gas of "compactons" is stable against 
arbitrary (square-integrablc) perturbations of the order 
parameter, hence it represents a local minimum of the 
free-energy landscape. This result is crucial to qualify 
"compactons" as the relevant effective degrees of free- 
dom responsible for sclf-glassiness of the complex fluid 
mixture. Therefore, we arrive to a very elegant and intu- 
itive picture of glassiness, as the nucleation of a "gas of 
compactons" , each of which corresponds to a local mini- 
mum of the free-energy associated with the LG eqs. |l]- 
[2]). Most remarkably, these "compactons" can be added 
together, each collection of "compactons" correspond- 
ing to a distinct dynamical partition of physical space. 
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This provides a very poignant and direct map between 
the complexity (coexistence and competition of multiple 
minima, sometimes referred to as " ruggedness" ) of the 
free-energy landscape and the morphological complexity 
of the fluid density in configuration space. This picture 
is highly reminiscent of the inherent-structures discov- 
ered/proposed based on numerical and analytical studies 
of glass- forming fluids [l(|. However, the link between 
the landscape and configurational complexity emerging 
from the present compacton-based LG picture, appears 
to be new and more direct. The same considerations ex- 
tend to higher dimensions, to be described in a future, 
more detailed publication. 

Since the collective properties of the "gas of com- 
pactons" shall be demonstrated via numerical simula- 
tions, it is crucial to prove that compactons survive dis- 
creteness, as we shall show in the sequel. In particular, 
we analyze under what condition on Do and D2 one can 
still find compact solutions on a lattice. To this aim, we 
considered stationary solutions of the discretized version 
of (fl} which become at x = and found a symme- 
try in the solution, namely x — > —x implies (j> ~4 > - 
This symmetry is clearly broken by the solution defined 
in ([5]) and the condition for the existence of a non-zero, 
symmetry-breaking solution of the discrete LG equation, 

reads Dp — + 2D2E > 0, Ax being the lattice spac- 



ing In the limit of small Ax and large D2, the 



latter yields D £ x -i < E and can be rephrased in terms 
of competing scales, as < 1, where is a scale 

proportional to Do/ Ax and Id has been defined previ- 
ously. This way, the limit Dq — > 0, where the system 
shows self-glassiness, reads as Id S> ^o- We now proceed 
to show that such self-glassiness is indeed observed in nu- 
merical simulations of the generalized LG eqs. ([T][2]). To 
this purpose, we simulated the generalized LG equation, 
including a noise term, to represent finite-temperature ef- 
fects. The corresponding Langcvin equation is simulated 
on a square lattice of size 256 2 with periodic bound- 
ary conditions. Initial conditions are chosen randomly, 
4>(x, y\ t — 0) = r, where r is a random number uniformly 
distributed in [—0.1,0.1]. In figurc[2l wc show two color 
plates of the order parameter <f>(x, y; t) at t = 20, 000 for 
the case Do = 0.3 and D2 = (top), D2 = 2.0 (bottom), 
both without noise. It is apparent how the case with 
D2 > leads to a much retarded coarsening, as a matter 
of fact to a dynamical arrest. 

In figure El we report the free-energy F(t) + 1/4 for 
Do = 0.6 and three different values D2 = 0, 4, 8 with 
e = 0. Each point is the result of the averaging on 
100 configurations with randomly chosen initial condi- 
tion. From this figure, it is seen that the asymptotic 
decay is always a power-law F(t) + 1/4 ~ t~ a , with 
an exponent, a, which becomes smaller and smaller as 
D2 is increased. Eventually, a(Z?2) reaches the zero- 
point (see the inset), formally corresponding to struc- 
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FIG. 2: Color plates of the order parameter <j>(x,y\t) for the 
case Do = 0.3, D 2 = (top) and D = 0.3, D 2 = 2 (bottom), 
e = 0, at t = 20, 000. The much slower coarsening associated 
with the D2 = 2 case is well visible. 



tural arrest, for D2 ~ 10.7. Next, we performed fur- 
ther simulations by including an external forcing, h, con- 
stant in space and time, as well as a thermal noise. 
We monitor the average response to the external drive, 
$(t) = M" 1 L- 2 YZ=i E x , v <t>m(x,y]t), where M is 
the number of realizations corresponding to different ran- 
dom initial conditions. With D2 = the system reaches 
its driven steady-state, <f> « 1 — 0(h) in a finite-time t c . 
As D2 > is switched on, this relaxation-time increases 
considerably Structural arrest similar to the one ob- 
served in figure El has been observed also in previous lat- 
tice Boltzmann simulations with full hydrodynamic 
interactions and conserved parameter dynamics. Figure[3] 
shows the relaxation time as a function of D2, t c (D2, Do), 
for Do = 0.3 and Do = 0.6 (inset) and e = 0.01. From 
this figure, it is seen that, as the ratio D2/ 'Dq is increased, 
the relaxation time starts to ramp-up quite rapidly. This 
divergence is consistent with a Vog el-Fulcher-Tammann 



law t c (D2,Do) = exp( D2 __ D2 ) [12|, where I?2,c and C 
both depend on Do and D2 plays the role of a temper- 
ature. In particular, we obtain Z?2.c ~ 2, and i?2,c ~ 12 
for Dq = 0.3 and Dq = 0.6, respectively. This ultra-low 
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FIG. 3: Free-energy decay for the case e = 0, Do = 0.6 and 
D2 = 0,4,8. The inset reports the exponent a(Z>2) of the 
corresponding power-law decay for D2 = 0, 2, 4, 6, 8. 
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FIG. 4: Divergence of the relaxation time t c at increasing 
values of D2, for Do = 0.3 and Do = 0.6 (inset). The noise 
amplitude is e — 0.01. 



relaxation is in line with the picture of a structural arrest 
of the mixture, due to the stability of the "compactons" . 
Another typical signature of glassy behaviour is ageing, 
i.e., the anomalous persistence in time of density-density 
correlations. A typical ageing indicator is the density- 
density correlator 



c(t w , 



((<l>(x,y,tw)<l>(x,y;t))c 
{4>{x,y;t w )^{x,y\t w )) t 



where t w is the waiting time and brackets denote spatial 
and ensemble averaging and (...) c stands for connected 
correlation. In figure [SJ we show this quantity for the 
case Dq = 0.6 and D2 = 4 and D2 = (inset) and 
e = 0.0003. From this figure, it is apparent that for 
D2 = the density-density correlator decays to zero, 
indicating that the system is able to visit all regions of 
phase-space. Such capability, however, is manifestly lost 
in the case D2 — 4, to an increasing extent as t w is made 
larger, which is precisely the ageing behaviour mentioned 
above. 

Summarizing, we have presented a new phase-field 
model exhibiting typical signatures of self-glassiness, such 
as long-time relaxation, ageing and long-term dynamical 
arrest. The distinctive feature of the present model is to 
allow the cost of building an interface to become locally 



FIG. 5: Time decay of the density-density correlator for var- 
ious values of the waiting time t w and e = 3 x 10~ 4 . The 
ageing effect, namely a decreasing loss of memory at increas- 
ing t m , is clearly visible. In the inset we report the case with 
D2 — for comparison. 



zero, while preserving global positivity of the overall sur- 
face tension. Analytical solutions are shown to take the 
form of compact density configurations ( "compactons" ) , 
associated with local minima of the corresponding free- 
energy functional. Direct simulations of the model show 
that self-glassiness emerges as a collective property of this 
"gas of compactons" . The compacton picture proposed 
in this work provides a very elegant and conceptually new 
link between the complexity of the free-energy landscape 
and the morphological complexity of the fluid density in 
configuration space. Valuable discussions with M. Cates, 
G. Gompper and I. Procaccia are kindly acknowledged. 
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